NMR properties of a one-dimensional Cu-O model 
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We obtain the Knight shifts and the relaxation rates related to the Fermi contact interaction term for 
a one-dimensional Cu-0 model using bosonization technique. We consider the small interaction limit 
at half-filling and away from half-filling. In this framework we predict that the antiferromagnetic 
contribution to the relaxation rate of the nuclear oxygen spin is completely suppressed even away 
from half-filling, when the temperature is low enough. In the strong interaction limit at half-filling 
we compute the effective Fermi contact interaction performing a Gutzwiller projection. Both limits 
suggest that the one-dimensional versions of the Mila-Rice and of the Shastry scenarios of transferred 
hyperfine couplings which were proposed to explain the NMR measurements for High-T c cuprates 
fail in a one-dimensional situation. 
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I. INTRODUCTION 

One of the challenging characteristics of the cuprate 
materials is the importance of magnetic fluctuations and 
their effect on normal state transport and on supercon- 
ductivity. This feature is naturally present in theoretical 
approaches emphasizing strong interactions in 2D. The 
role of magnetism can also be assessed in other scenarios 
promoting marginal or nearly antiferromagnetic Fermi 
liquid (NAFL) behavior. The N AFpL framework has been 
used by Millis, Monien, and Pineail to diacuss NMR e: 
periments: on the basis of the Mila-Ricen and Shastr 
local terms describing spin fluctuations induced by the 
hyperfine interactions, these authors were able to com- 
pute various Knight shifts and nuclear relaxation times. 
However, there are critical remarks and further details to 
this theoryoQ and also other theories introduced in or- 
der to interpret the behavior of the various Knight shifts 
and nuclear relaxation times in high-T c superconductors, 
which we do not discuss hereB 

In this paper, we will follow the basic assumption of 
Mila and Rice that it is necessary to include a sizeable 
isotropic hyperfine interaction term to fit the data of 
NMR experiments in High-T c cuprates. Thus in the 
following we will focus our interest on this contribution 
which comes mainly from the Fermi contact interaction 
between a nucleus and its surrounding partially filled s 
orbitals, namely the 4s orbital for copper and the 3s or- 
bital for oxygen. Experiments have shown that the mag- 
netic properties are described by a single-spin component 
model. In the weak interaction limit, this single-spin de- 
gree of freedom could be associated with tb&. strongly 
hybridized Cu-3d — 0-2p anti-bonding band,B0 whereas 
in the strong interaction limit, it is associated with the 
nearly localized Cu-3d spin.l3'EH3 Within the local picture, 
the contrasting NMR behavior seen on the Cu and O sites 
arises from their different hyperfine form factors. These 
are nothing but the Fourier transforms of the Fermi con- 
tact interaction terms approximated by a sum over sur- 
rounding localized next-neighbor Cu-3d spins. 

One of the (many) complications concerning the 



physics of cuprate materials is that there is still no con- 
sensus about what should be the correct theoretical ap- 
proach to treat correlations in 2D: Is the ground state 
Fermi liquid or non-Fermi liquid like? Can one treat 
interactions perturbatively, or is it more appropriate to 
treat kinetic terms as corrections in the strongly interact- 
ing limit? In the latter category, working out a consistent 
treatment of the non-double occupancy constraint is still 
an open issue. 

By contrast, one-dimensional systems offer a perfect 
testing ground for the study of magnetic fluctuations. 
Since in one dimension it is possible to treat correlation 
effects properly both in the limit of weak and strong in- 
teractions, such models allow to compute explicitly the 
dependence of the relaxation. This allows to get some 
feeling for the effects of doping. In addition to the insight 
that such study allows to gain for higher dimensional 
models, there are explicit realizations of one-dimensional 
systems, suck-as the Bechgaard salts or copper germanate 
compounds .Bllj 

For these reasons, we choose to investigate hyperfine 
interactions in the one-dimensional version of the Cu-0 
model. This allows us to extract form factors both in the 
insulating and in the doped regime, without assuming a 
specific form of the Fermi contact interaction term. We 
can then compare the exact results with the predictions 
that the standard approximation schemes used in two 
dimensions would give in the one-dimensional situation. 

The paper is organized as follows. In section ||, we in- 
troduce the model in one dimension, as well as the three 
different approximations used to describe the magnetic 
relaxation processes related to the Fermi contact interact 
tion term , namely that due to Mila-Rice,B to ShastryEl 
and to Bulut.Q In section III, we solve the full model for 
weak interactions as compared to the bandwidth. We 
obtain spin-spin correlation functions and discuss the 
Knight shift K and the relaxation rate 1/Ti in detail. 
We compare our results with the prediction of the Bulut 
model, which is applicable for weak interactions. Sec- 
tion [Fv| solves the problem in the opposite limit of very 
strong interactions using a Gutzwiller projection climi- 



nating double occupation on the copper sites. We again 
compare our results with the one-dimensional extrapo- 
lation of the Mila-Rice and of the Shastry approxima- 
tions. A general discussion of our results as. presented 
in section M. Since NMR data on organic^ and inor- 
ganic quasi-lD compounds seem to give an essentially 
isotropic relaxation rate, the body of the paper mostly 
focuses on the isotropic contribution to the hyperfine in- 
teraction. Yet, for the sake of completeness and in view 
of the fact that K can be anisotropic (see below in section 
V) we discuss the effect of anisotropic hyperfine terms in 
Appendix H these terms only modify prcfactors in the 
expressions of \jT\ and of K. Appendix M and offer 
details of our calculations. 



II. DEFINITION OF THE MODELS 

A. The four-band model 

We consider a system with two different atoms per unit 
cell (Cu and O). In order to describe the ground state 
properties, we take into account the 3d and 2p orbitals 
on Cu and O respectively. The related hole states are 
represented in Fig. tt\ and denoted by a and b in the fol- 
lowing. Since the coupling to the nuclear spin via the 
Fermi contact interaction occurs only for partially filled 
s orbitals, we also have to retain the Cu-4s and 0-3s 
shells (denoted by A and B) to correctly obtain the de- 
sired NMR properties. The Hamiltonian describing the 
system can thus be written as 



H — Hn + H,Q 



H 



N 



(2.1) 



where H contains the electronically relevant orbitals a 
and b. H$ describes the coupling of orbitals a and b to 
orbitals A and B. As indicated in Fig. [y, the orbitals A 
and B are basically filled and produce only small correc- 
tions to the electronic term represented by Hq , so that we 
will treat H$ as a perturbation. Finally, Hpj describes 
the coupling of the orbitals A and B to the nuclear spins 
and will be treated as a small - time dependent - per- 
turbation in linear response. 

The main contribution, H$, is given by 



Hq = Ht + Hi 



(2.2) 



where 



Ht = J2 e - 7 



e b nbj 



~ y^afe 0/ ja .(bj a + bj-l >a ) + h.C. 
ja 

Hu = ^ U a n a ^n aji + ^2 Ubn bj ^n bj 



(2.3) 
(2.4) 



Here t a b describes the hopping between the Cu-3d and 
the 0-2p orbitals with the phase conventions shown in 



Fig. g. U a and Ub are the local repulsions on the cop- 
per and oxygen sites as shown in Fig. |lj. The Coulomb 
repulsions Ua and Ub can be ignored assuming that en- 
ergy cost considerations discourage processes in which 
two holes are excited in one s orbital. A nearest neigh- 
bor interaction U a b could also be added to the model to 
generate a phase transition to a superconducting phase,E3 
but here our study deals with the vicinity of the half-filled 
case, i.e. near the antiferromagnetic phase, and we will 
ignore U ah . 

The coupling between the orbitals a, b and A, B reads 

H S = ^ e A n A + ZBUBj 

j 
+ 22 t Ba Bj a (aj+i t<T - a ja ) + h.c 

jo- 

+ J2 t Ab A]„{b ja + bj-i t( r) + h.c. 
+ 5Z tAa A' jcr (aj + t, a + a,j-i >a ) + h.c 



(2.5) 



jo- 



with the pha se c on vent ions o f Fig . [2| The density oper- 
ators n v j in (|2.3|), (|2.4p, and (pTq) are the standard ones 



(2.6) 



/ , n r\ja — 2_^i T ljcr r ljo 



where r\ = a, b, A, B. 

Finally, the isotropic coupling to the nuclear spins (I 
for the copper atom and J for the oxygen atom) is given 
by 



Hn — 2_^ CaIjSaj + CbJjSbj , 



(2.7) 



where similarly to (2.E) the spin operators are given by 
S '"' = 2 51 4<yi a o 1 a 2 ri ] a 2 ■ (2.8) 

(T\IJ2 

The coupling constant for the electron-nuclear interac- 
tion is given by C v = (8Tr/3)\ip ri {0)\ 2 jn9^B^, which is 
proportional to the local hole density for the respective s 
orbitals at the origin. a ai(72 are the Pauli matrices. 

B. Reduced models 



In ID, the Hamiltonian (|2.l| ) can explicitly be ex- 
pressed in terms of bose fields. This allows for a full 
treatment, in which Hq, Hs, and Hn are all treated on 



an equal footing. This is the route we follow in section III 



In higher dimension there is still no solution of the fully 
interacting problem. So various approximation schemes 
have been devised and applied to each of the pieces of 
H separately. These lead to effective " nuclear" Hamilto- 
nians : the Mila-Rice and the Shastry model for strong 



interactions and the Bulut model for weak interactions. 
Typically, the analysis by Mila-Rice and Shastry starts 
from a partially projected Hamiltonian 



H = PH Q P + H s 



H 



N 



(2.9) 



The first part gives the t-J model or at half-filling the 
Heisenberg model, which contains the dynamics related 
to the Cu-3d and 0-2p orbitals. The second part con- 
tains the unprojected degrees of freedom related to the 
Cu-4s and 0-3s orbitals, as well as the electron-nuclear 
interaction part. Further approximations for Hg and H^ 
lead to the Mila-Rice model or to the Shastry model (see 
below). In the weak interaction limit, Bulut et al. have 
proposed an RPA treatment of Hq in combination with 
the effective electron-nuclear interaction term of Mila- 
Rice and Shastry. 

Before we turn to the full solution of the model, let us 
review the main features of such approximations when 
applied to our one-dimensional system. This will allow 
us to contrast the predictions of the ID version of these 
three models and the results obtained for the four-band 
model, which may provide some clue to the validity of 
these approaches for strongly correlated systems. 



1. The Mila-Rice model for strong interaction 



The model defined in (2.1) is approximated by 



H~Hr + H N 



Mi 



(2.10) 



Hq % is the approximation for ( |2.2| ) and denotes, at half- 
filling, a Heisenberg model for local Cu-3d spins gener- 
ated by 



ff Mi = P[H (U b = 0)]P , 



(2.11) 



where P is the Gutzwiller projection operator which pro- 
hibits doubly o ccup ied Cu-3d states. The additional un- 
projected part ( j2.5| ) with tAa = ts a = an( i the electron- 
nuclear interaction part (|2.7|) are approximated by 



TjMi 



Mi 



^2C A ljS Aj 



where 






F Aa ( S aj"-1 + S 



o,i+i/ 



(2.12) 



(2.13) 



Thus S A il is the Mila-Rice approximation for the original 
spin Saj used to explain the NMR experiments measured 
on the copper sites. F%£ = \X A ^\ 2 denotes the effective 
overlap between one Cu-4s spin with a neighboring Cu- 
3d spin. In the Mila-Rice model only hopping processes 
via the 0-2p orbitals are included, whereas the direct 
hopping between Cu-3d and Cu-4s orbitals is ignored. 
In Ref. @ Mila and Rice perform a quantum chemical 



analysis without including interaction effects, so we will 
do the same and propose for the amplitudes 



\ Mi 
A Aa 



tAatab 



(e a -e A )(e a - e b ) 



(2.14) 



This result is obtained by a projection in real space of a 
Cu-4s orbital onto a neighbo ring C u-3d orbital for U a — 
0. The Fourier transform of (2.13) is given by 



(2.15) 



(2.16) 



&Ap — F Aa(.P)^ap 

with the Mila-Rice form factor 

F^(p) = 2F^cos(pa). 

2. The Shastry model for strong interaction 



The approximation proposed by Shastry in Ref. |3| is 
given by 



where 



U ^ uSh I rrSh 



H$ h = P[H (U b = 0,U a = oo)]P 



(2.17) 



(2.18) 



leads to the Heisenberg model at half-filling and to the 
t-J model for a doped system with strong repulsion on 
the copper sites. The electron-nuclear interaction part 
reads 

H N = } y CAljS Ai + C B JjS B j , (2-19) 



where the spins S^' 1 are approximated by a linear com- 



bination of unprojected Cu-3d orbitals 

iSh 
•>Aj 



SSh rpSh fc i C ^ 

Aj — t Aa (Pa,j-1 + Ofl,j+lJ 



+ F Ab (Sftj +Sb,j-l) 




— F Aa ( S a,j-1 + S a,j + l) 


(2.20) 


and 




& Bj — ^Ba (Z>aj + »aj + lj 


(2.21) 


with the coefficients 




F S A h a =\Xf a ?=( ^ V 
\e A - e a J 


(2.22) 


f% = \\aI\ 2 = (^-) 2 

V e A — £fc / 


(2.23) 


TpSh \\Sh\2 ( tBa \ 
b Ba = lABal = , 


(2.24) 



es - (a. 



For finite doping, it is assumed in Ref. that the spin de- 
grees of freedom related to S^ are quenched in a Zhang- 
Rice singlet.E3 This a ssump tion justifies the second ap- 
proximation done in (2.20). Further, Shastry includes 



only the direct couplings ( 2.22J ) up to second order pro- 
portional to t Aa for the relaxation of the nuclear copper 
spin and ignores the fourth order contributions propor- 
tional to t 2 ab t 2 Ab via the 0-2p orbital as proposed by Mila 
and Rice. The Fourier transform of the approximated 
spins S:? h reads 



cSh o j^Sh 



with 



2Op)S 



F^(p) = 2F^cos(pa) 
F^(p) = 2F§ h a cos(pa/2) 



(2.25) 

(2.26) 

(2.27) 



For the uniform contribution (p ~ 0) all form factors 
are finite, but for the antiferromagnetic wave vector 
(p ~ 7r/a) the form factor vanishes for the oxygen sites, 
where as it s tays finite for the copper sites. Thus (2.16), 



( 2.26 ), and ( 2.27 ) are the one-dimensional analogs of the 
NMR Mila-Rice and Shastry form factors for High-T c 
cuprates. 



III. THE WEAK INTERACTION LIMIT 

Let us now solve the full model ([2^) when interactions 
are weak compared to the bandwidth. This allows us to 
use the bosonization technique for treating interactions 
in the undoped as well as in the doped case. 

A. NMR properties of the four-band model 

1. Reduction to an effective single-band Hamiltonian 

Instead of working with the basis a, b it is more con- 
venient to diagonalize (|j|) within a unit cell, and to 
introduce the bonding and anti-bonding bands. Using 
the transformation 



aka = [cos(7 fc )a fe(T - sin(j k )f3 ka ] e 
b ka = sin(j k )a ka + cos(j k )/3 krT 



(3.1) 



with 



3. The Bulut model for weak interaction 

In Ref. 0, Bulut et al. used a weak interaction RPA 
calculation combined with an electron-nuclear interaction 
as proposed by Mila-Rice and Shastry to compute NMR 
related quantities for the High-T c cuprates. The one- 
dimensional analog reads 



H~H, 



Bu 



ttBu 



(2.28) 



In two dimensions, H$ u is obtained by applying the RPA 



method to the 2D version of (2.2) with XJ\, = 0. In one di- 
mension, we can treat all the interaction terms (U a ,Ub) 
of the original Hamiltonian Hq by means of bosoniza- 
tion and of renormalization-group theory. Finally, the 
electron - nuc lear interaction term H^ u is given in anal- 
ogy to ( [2.19J ) replacing the approximated spins by 

Fa: (S«j-i + Sa, j+ i) (2.29) 

F§Z(S aj + S a , j+1 ) . (2.30) 

In Ref. p], the parameters F^" are undefined and could 
in general include all possible overlaps of the Cu-4s and 
0-3s orbitals with the Cu-3d orbitals in the sense of Mila- 
Rice and Shastry. The coefficient for the oxygen will be 



o5u 
a Aj 

cBu 
»Bj 



F, 



Bu 



|A 



Bu\2 
Ba\ 



\\Sh\2 
\*Ba\ 



(2.31) 



and is thus the same as that proposed by Shastry, 
whereas the coefficient for copper 



* Aa 



I \Bu\2 
\*Aa\ 



|A 



Sh | \Mi\ 
+ A Aa\ 



(2.32) 



includes additional combined terms of third order pro- 
portional to tAatAbtab, which are absent in the Mila-Rice 
and the Shastry model. The form factors correspond to 



F£(p) = 2F£coB(pa) 
F^(p) = 2F^co S (pa/2) 



(2.33) 

(2.34) 



tan(27 fe ) = —2- cos(ka/2) , j k e 
e 

the kinetic energy ( |2.3| ) becomes 



ha 



(3.2) 



(3.3) 



where the state \a ka ) refers to the lower Hubbard band 
with energy e a (k) — —e/cos{2'y k ), and the state \0ko) 
to the upper one with energy ep(k) = e/ cos(27fe). In the 
absence of interactions the chemical potential /i lies in the 
a-band both for the undoped and for the doped system, 
and one can ignore the /3-band, which is at least tb — e a = 
2e higher in energ y. T he same prope rty holds when the 
interaction terms (2.4) are added to (EOj), given that in 
the weak coupling limit Ub,U a <C 2i^/e. Correlation 
effects in ID will strongly affect the a-band states, thus in 
the following we ignore the terms containing /3-operators 
when substituti ng ( | 3.1 |) into (2.4). 

Substituting ( |3.1[ ) into (2.5) and performing a first or- 
der perturbation theory with respect to H$, all operators 
in (|2.1|) can be written as 



with 



Vka — \ a (k)a ko 



X aa (k) = cos(7 fe )e~ 1 ^ 

X ba (k) = sin(7 fc ) 

2i^ a cos(fca) . 
cos(7fc)e 



(3.4) 



XAa(k) = 



e a {k) - e A 



2tAbCOs(ka/2) . , _,ko 

+ m sm( 7fe )e * 

e a (k) - e A 

, , n 2it Ba sm(ka/2) 

Abq(k) = — cos(7fc) . 



(3.5) 



Here we have assu med an unperturbe d g round state 



a (1 >) = \a). Thus (|3.4| ) implies that (|2.l|) reduces to 



an effective single-band Hamiltonian. 



2. The Continuum limit 

We can now use the standard techniques in order to 
treat interacting one-dimensional systems. Restricting 
ourselves to the low energy physics regime we make the 
usual approximation valid for ID systems, i.e we lin- 
earize the spectrum close to the Ferm i points, as shown 
in Fig. g. Then the Hamiltonian (2J2) is reduced to 



H T = ^Z rv F qal qcr a rqa (3.6) 

r=±,q,a 

Hu = Yl -^ a ti, q i+m1 a t2, q 2- q 3l a ^^iar 4 , qi ^ , (3.7) 



r.q 



where r = (ri,r 2 ,r 3 , r 4 ) and q = (01,02,03)- U(r) pa- 
rameterizes the repulsive interaction in the continuum 
limit and is given in terms of the standard notations as 

U(±,±,±,±) = U 

U(±,T,±,T) = U 8 

U(±,±,T,T) = U C . (3.8) 

U refers to the two forward scattering processes, U s to 
the backward scattering, and U c to the umklapp scatter- 
ing process that occurs at h alf-fi lling. The relation to the 
local repulsions defined in (2.4) is given by 



Uo = U ba + U aa >0 
u s = U ba + U aa >0 
u c = U ba -U aa <0 



(3.9) 



where V^a is the short notation for the projected 
Coulomb energies U r} \X r]a (kF)\ 4 ■ 

As usual for interacting one-dimensional systems, it is 
useful to introduce a boson representation of the fermion 
operators, related to the charge and spin density fluctua- 
tions. Since the technique is standard, we only r£caLL the 
main steps and refer the reader to the literatureolHI We 
rewrite the original density operators in terms of a linear 
combination of charge (y = c) and spin {y = s) density 
operators for each branch 



Pr 



(Pr 



ap, 



0/V2- 



(3.10) 



These density operators define the phase fields 
in v-^ 1 



Ms) 



L ^„ q 



r,qjt0 



ITT v— v r 



e„(*) = y £ 



a\q\/2~iqx 

O fJ rL / 



-a\q\/2 — iqx 



L " q 



(3.11) 



All op erators can be expressed in terms of the boson fields 
( |3.11 ), and the fermion operator reads: 



a ra (x) = 



1 



/ 27ra 



ivkpx 7= 

=e V2 



>(* c +<7* s )-(e c +<Te s )] 



The complete Hamiltonian becomes 
where 



m= ^ 



dx 

2^ 



{u v K v ){-KU v f + [%- W,) 



K, 



(3.12) 



(3.13) 



(3.14) 



is a quadratic part containing o nly ch arge or spin degrees 
of freedom (with v = c, s). In (3.14), the variable H^ = 
d x Q v is the momentum density conjugate to <&„, and thus 
they respect the commutation relation [^ ll (x),JI v (x')] — 
iS(x — x'). The interaction terms are given by 



m 



m 



dx— ^r cosfVs^c — Sx] 



dx 



(2™) 2 
(2ira) 2 



cos[V8$ c - S(x - a/2)] 



, 2aU s , f-* 
dx (2^ C0S[ ^ $s 



(3.15) 
(3.16) 



Here 8 = 4kp — 2n/a is proportional to the doping of 
the system with respect to the half-filled case shown in 
Fig. p] (for which Uf = it /2a) . Using this representation 
we suppose to work with a fixed number of particles, since 
kp is directly related to the filling. Finally, the isotropic 
electron-nuclear interaction part could be written as 

H N = f dx[C A I(x)S A (x) + C B J(x)S B {x)} , (3.17) 



wher e Hn is th e projection of (2.7) onto the a-band using 
( |2.8|) and (3^). Th e pro jected spin operators S^ are ex- 
pressed in terms of ( 3.1 2| ) . For example, the ^-component 
of the spin operator S,, can be represented as a sum of 
p ~ and p ~ 2k p components as 



S%(x) = |A^ Q (fc_F)| 2 [s Q (x) +s va {x)] 
where the non-oscillatory part is given by 

1 



s a (x) 



-{d x $ s ) 



V2tt 
and the oscillatory part by 

s V a{x) — — sin[v2$ s ] sin[2fc_F(a 



(3.18) 



(3.19) 



V2$ c 



(3.20) 



The difference between the copper and the oxygen sites 
is reflected in the value of x n and affects the oscillatory 
part; indeed, for copper x a — xa = a/2 and for oxygen 



Xb = %b = as a consequence of the different phase 
factor s in ( 3.5 ) . 

In (3.14), the u v are the new velocities for the v- 
excitation and the K v are the Luttinger liquid param- 
eters controlling the anomalous exponents in the corre- 
lation functions. F or w eak coupling, they are related to 
the interactions in (|3.7j) by 



u s K s = u c K c = Vf 
u s /K s = Vf - aU /n 
u c /K c = vf + aUo/n 



(3.21) 



Since the Luttinger liquid representation is more general 
than the perturbative result for small interactions, it is 
also applicable when the interactions are strong. The 
quadratic Hamiltonian can be viewed in this case as an 
effective Hamiltonian describing the low-energy proper- 
ties of the system, provided that the correct Luttinger 
liquid parameters are used. Such a smooth connection 
between weak and siipng coupling has been proven for 
single-band models JSEj and a similar Luttinger repre- 
sentation has-bpen shown to work f or th e case of the two- 
band model.EHlEj Equations ( |3. 13 -3. 17) define the four- 
band model, and the NMR properties can be computed 
through Hn. 



3. Correlation functions at zero temperature 

We focus here on the spin-spin correlation functions 
relevant for NMR and for neutron scattering experi- 
ments. The general form of these functions is 

R w (x,t) = (T t S*,(x,t)S*(0,0)) - (S z n ,) (S z v ) , (3.22) 

and it describes correlations between different orbitals rj 
and rj' at different points in Euclidean space-time. Here 
we introduce the decomposition of this function into a 
non-oscillatory and an oscillatory part 

linn' (x, t) = R, m > (x, t) + cos(2k F x)R r)r)l {x, r) , (3.23) 

since the behavior of these functions will be very different 
for one-dimensional systems. 

Because of the doping dependence in the cosine terms 
in ( |3.15| ), the behavior of the system will quite clearly 
be different for zero and for finite doping. At half-filling, 
one sees from ( |3.21 ), (B2), and (B3) that charge exci- 
tations are massive (c m ), whereas spin excitations are 
in the massless regime (s ). One recovers the standard 
Mott or charge-transfer insulator with the massless exci- 
tations corresponding to a Heisenberg-like exchange. In 
the doped case, the term ( 3.15J ) is irrelevant because of 
the oscillatory factor Sx. However, at short distances 
or for short times this term is still small, and the co- 
sine term will influence the behavior of the system. We 
thus distinguish between two different regimes for the 
doped case: we assume that for intermediate distances 



(a « i « Is) the system remains in the (c m , s )-phase 
as mentioned before for the half-filled case, and when dis- 
tances are larger than Is, the system will be in the (c , s )- 
phase because the umklapp process becomes ineffective. 
The characteristic length separating these two regimes 
denotes essentially the distance between two charge do- 
main walls and is given by Is = 2ir/5. 

Due to the spin-charge separation in ( p,13| ), each part 
of the correlation function ( 3.23| ) will factorize into in- 
dependent averages over the spin (s ) and the charge 
sector (c D or c m ) , and will only depend on the charac- 



teristic distance r v 



[(UvTY 



between two points 



in Euclidean space-time (with v = s,c). Details about 
the correlation functions in the various regi mes (cj, Sj) 
are expl ained in Appendix 0. Substituting (3.18-3.20) 
in ( 3.22 ), the non-oscillatory contribution to the correla- 
tion function is given by 



^i]rj' I ^ 



<K» 



j |At7'q,| R a \Ts) 



(3.24) 



where R a {r s ) — {2^r s )~ 2 depends only on the spin de- 
grees of freedom and is thus completely independent of 
the coexisting charge phase. Notice that this function is 
also independent of the orbitals rj and rj' , and thus there 
is no fundamental difference between copper and oxygen 
contributions. 

For the oscillatory part of the spin-spin correlation 
functions, the situation will be quite different. We re- 
strict ourselves to the calculation of correlation functions 
between identical orbitals (rj = rj'). Using average s ove r 
the charge and spin sectors of the Hamiltonian ( 3.13 ), 
these functions can be reexpressed as 



jbc m ,s 



\\a\ ' 



R%(r c )R s a °(r s ) 



(27ra) 2 
in the massive charge regime and as 



hc 0: s _ I^Jjal 



(2na)< 



R c a °(r c )R s a °(r s ) 



(3.25) 



(3.26) 



in the massless charge regime. The newly defined cor- 
relation functions in the massless phases (z/ ) are given 
by 



K°{r v ) = (a/r u )"»F(r v ) . 



(3.27) 



The function F(r„) describes the corrections to the Lut- 
tinger Liquid behavior which come from the flow to the 
fixed pointEil. To lowest order, F{r,j) can be approxi- 
mated by 1. The renormalized Luttinger liquid parame- 
ters K* for a spin symmetric model with repulsive inter- 
action arc restricted to 



K* = 1 and < K* < 1 



(3.28) 



The value of the renormalized Luttinger liquid parameter 
K* depends on the interactions. For weak interaction, 
K* is close to 1, and it decreases as interactions become 
more repulsive. 



The correlation functions in ( |3.25| ) which are charac- 
terized by the massive charge phase are given by 



R< A a ( r c) = 2cosh[K c K (m c r c )}(m c a) K " 
R C B a a ( r c) = 2smh[K c K (m c r c )](m c a) K ' : 



(3.29) 



5. Knight shifts and relaxation rates 

The standard expressions for the Knight shifts and for 
the relaxation rates resulting from a hyperfine coupling 
term like (2.7) are 



and depend on the chosen orbital 77. Thus, the behav- 
ior for copper and oxygen will be quite different. It de- 
pends on the distance r c , the mass m c , and the stiff- 
ness constant K c . In general, for distances larger than 
lm e = V TO c, the function R C J^ for copper tends to a finite 
constant, whereas R C g l a for oxygen tends exponentially to 
zero. 



K 1 



(->, 



1 

t7 






It? Jr,left (3 ^ 



2^ 



(3.33) 

(3.34) 



where 7 = refers to the half-filled case and 7 = 8 to the 
doped case, uj„ denotes the electronic Zeeman frequency 
in orbital 77, which is very small as compared to the en- 
ergy scale of the purely electronic system fixed by the 
cutoff A. For the Knight shifts the sum is restricted to the 
active orbitals a and b. We can split up the susceptibil- 
lation function i^fejn/3), we will only use the asymp- ity % 7 into the non -oscillatory x 7 , and the oscillatory 



4- The asymptotic expressions at finite temperature 



In order to obtain the temperature dependent corre- 



totic expressions of ( 3.29J ) . We recover a Luttinger liquid 
behavior, and the temperature dependence can easUv be 
obtained with the help of the conformal symmetry;L3 in- 
deed, we only need to replace r v (x, r) by r v (x, r, 13) where 



r v (x,r,/3) 



u u (3 



' sinh 



X — lUyT 



Iv/3/tt 



sinh 



x + iu v T 



UvP/tT 



(3.30) 



The relevant asymptotic expressions at half-filling and 
away from half-filling depend on the relative magni- 
tudes of the various characteristic lengths of the sys- 
tem, namely the lengths related to the mass, l mc , and 
to the doping, Is, as well as the thermal length lp — 
mm{(u c /3) -1 ,(ti a /3) -1 }. For half-filling (7 = 0) at low 
temperature, we are in the regime where lp ~^> r c ~^> l mc 
and Is = 00, thus we can approximate the oscillatory 
charge contribution in (3.25), and the functions are sim- 
plified to 



R r,ri = C 7]W<*\ R a (r s ) , 



(3.31) 



where the amplitude of the oscillatory part at half-filling 
Rfi™ (oo c ). The remaining correlation 



is given by C° 



function is independent of the orbital 77 and given by 
i?° = (2na)- 2 (a/r s ). 

The large distance limit of the corresponding expres- 
sion for small doping (7 = 8) and low temperature, where 
^ > r c > Z5 > lm c , looks like 



Ri 



CAXrja] R a (r c ,r s j 



(3.32) 



In the doped regime, C^ is the amplitude R^(h) ob- 
tained in the massive phase at the crossover, as shown 
in Fig. [|. Like before, the remaining correlation function 
R 5 a = (2na)~ 2 (a/r s )(a/r c ) K <= is also independent of r\ but 
shows dependence on spin and charge degrees of freedom. 
For larger doping rates (Is < l mc ), the difference between 
copper and oxygen sites vanishes. 



contribution xl , just like for the correlation functions 



in ( |3.23 ). Finally, the Knight shifts for the linearized 
four-band model in units of C r/ /(7 r? 7 e fi. ) are given by 



K^=F r , a Xa(u> = 0,q^0) 



(3.35) 



and the relaxation rates in units of C 2 / (^f^Ji oj v ) by 



TpT = ~s E Im ( F n<*) 2 Xa(wr,,q) + (F^ a ) 2 xl(ur h q) 



T 



1'/ 



B 



\q\<\ 



(3.36) 



The susceptibilities Xa an d xL m space-time can be ob- 
tained from 

X a (x,t) = 20(t)Im[R a (x,r.j3)] T=it+e (3.37) 

Xl(x,t) = 29(t)lm[RZ(x,r,f3)] T=it+e (3.38) 

performing the continuation to real time. The r-ordered 
temperature dependent Green's functions on the right- 
hand s ide a re the same as in (3.24), (3.31) and ( 3.32|) , 
using ( 3.30 ). Thus, in general we can calculate ( 3.35 ) 
and fl3.36|) by performing the Fourier transform of (3.37) 
and ( |3.38| ). Here we restrict ourselves to the solutions 
obtained by the so called power counting method. The 
temperature dependences of the Knight shifts and of the 
relaxation rates are shown in Table H and the form factors 
F va and F^ a are given in Table ||. 



B. NMR properties of the Bulut model 

In order to obtain the NMR properties of the ID ver- 
sion of the Bulut model we perform the same procedure 
as before for the four-band model. The bosonized version 



of the Bulut model (2.28) is given by (3.13) replacing Hjy 
by 



H Bu = fdx [C A l{x)S B A u {x) + C B J(x)S% u (x)] 



(3.39) 



where S^ " is the proje ction onto the a-band using ( 2.29 ) , 
( [2.30 ), (2.£) and (3.4). For the z-component of the spin 



S% u we get 



S A Bu (x) 
S B Bu (x) 






2 [s a (x) + cos(2k F a)s aa {x)] 
2 [s a (x) + cos(k F a)s aa (x)] 



(3.40) 



ignoring all gradient terms of the field <j> c . Th e spi n op- 
erators s a and s va are defined as before in ( 3.19 ) and 
(1 3.20 ). The NMR p roper ties fo r a h yperfine coupling like 
( |3.39 ) are given by ( 3.35| ) and ( 3.36J ) by the replacements 



F, 



rja 



F, 



The values for the differ- 
ent form factors are shown in Table [II. At this level of 



1 1 a. 



and F^ a 



P7« 



approximation, both the four-band and the Bulut model 
show exactly the same temperature dependence for the 
Knight shifts as well as for the relaxation rates; this de- 
pendence is different for the uniform contribution and 
for the oscillatory one (see Table |), as is well known 
for interacting one-dimensional systems. This effect has 
nothing to do with the various orbitals where the Knight 
shifts and the relaxation rates are measured. 



C. Comparing the four-band model and the ID 
Bulut model 

First we focus on the coefficients of the Bulut and of 
the four-band model (compare Table || and III ) related 
to the different projection procedures of the s-orbitals 
onto the ground state. For comparing both models, we 
investigate the limit t a b <C (e& — e ). Then, for the four- 
band model the projection of the s orbitals (A, B) onto 
the lowest band (a) is strictly done in fc-space and results 



\\B a {k F )f 



2(Afj 2 
2cos(fc F a)(A^ a ) 2 



(3.41) 
(3.42) 



|A^(M| 2 ^2(Af a + A^) 2 (3.43) 

+ 4(A!*) 2 (3.44) 

+ 2cos(2fc F a)(A^ + A^) 2 (3.45) 

+ 8 cos(k F a)(Xf a + \%)\% , (3.46) 

whereas for the Bulut model it is a combination of real 
space and fc-space projection yielding 



2|A^| 2 



Sh\2 



|A OQ (Mr->2(A^J 



iiWi\2 



2|A^|^|A QQ (fc F )r-2(A^ a + A^) 



(3.47) 



(3.48) 



The general solution f or th e projected 0-3s orbital in- 
cludes one more term (3.42) than the solution proposed 



by Bulut (3.47). This term corresponds to a dynamic 
contribution which includes a charge displacement. How- 
ever, for a half-filled system the additional term van- 
ishes and the two solutions become identical. By con- 
trast, the projection procedure for the Cu-4s orbital pro- 
duces a completely different behavior in the two mod- 
els. For a half-filled system, the hopping processes via 
^Aa + ^Aa contribute only in the Bulut model 



\Bu 



(I 



48), whereas they are exactly canceled by the related 
dynamic terms ( 3.45 ) in the four-band model. Thus, for 
the four-ba nd m odel at half-filling, only an additional 
local term (3.44) as well as a dynamic combined term 
( |3.46 ) remain. The term (3.44) is the local analog to the 
t rans ferred terms proposed by Mila-Rice, and the term 
( |3.46|) is a combination of Mila-Rice and Shastry terms 
which includes a charge displacement. It should be clear 
that our projection procedure is the right one for a sys- 
tem with small Coulomb interactions: First we diagonal- 
ize the tight-binding Hamiltonian dealing with extended 
wave functions, and then we treat the Coulomb energy 
approximately within this non-local basis. The approx- 
imation proposed by Bulut suffers from a mismatch be- 
tween the local and the non-local point of view. 

The second part of the oscillatory contribution to the 
form factors (compare Table O and [II), which contains 
the dependence on the characteristic lengths related to 
the doping rate, Is, as well as to the charge mass, l mr , 
is the crucial one. Away from half-filling, the four-band 
model shows a different behavior on the copper and on 
the oxygen, despite the fact that umklapp processes only 
contribute on short or intermediate scales. Indeed, the 
different hyperbolic dependencies of the two character- 
istic lengths l mc and Is for copper and for oxygen (see 
Table ||) affect measured quantities related to long dis- 
tance or long time behavior. Instead, for the Bulut model 
the difference between copper and oxygen comes in only 
because of the special choice of a Mila-R ice-Shastry type 
electron- nuclear interaction term (3.39) and the related 
unconventional projection procedure which res ults in the 
different trigonometric form factors (see Table III). The 
influence of the charge mass m c is the same for copper 
and for oxygen, a fact which manifests itself by the same 
dependence on the characteristic length l 7ric . 

Note that the four-band model leads to a very small 
contribution on the oxygen even at finite doping, be- 
cause the contribution is exponentially suppressed in a 
way which depends on the ratio between l mc and Is, 
whereas the oscillatory contribution on the copper atom 
is nearly independent of the doping rate for long distances 
or times. For the ratio between copper and oxygen we 
distinguish between two regimes: 



tanh 



K C K 



for 
for 



IS < It, 

IS > l r . 



(3.49) 



In the former regime, we recover the Luttinger Liquid 
behavior, since the infinite length l„ la stems from the 



vanishing of the umklapp process when U ac 



U bo 



that case there is no fundamental difference between cop- 
per and oxygen anymore. Only the overlaps with the 
ground state remain different. The latter regime, where 
the fundamental difference occurs, will be reached ex- 
ponentially as K c y^7rl m< ,/2ls exp(— ls/l mc ), and thus the 
oxygen does not see the antiferromagnetic fluctuations 
in this limit. Instead, for the Bulut model everything 
depends on the same correlation function, and the dif- 
ference between copper and oxygen comes from the fil- 
tering factors. Thus, the oscillating contributions to the 
relaxation rates for the oxygen is always proportional to 
[0 + (Tra/2ls) 2 ], whereas the contributions on copper are 
reduced by a factor [1 — (na/ls) 2 }. The ratio of the os- 
cillating contribution to the relaxation rates is approxi- 
mately given by (ira/2l$) 2 , and is completely independent 
of the details of the projected local Coulomb repulsions 
U aa and Ub a - It only depends on the doping rate and 
is proportional to (Sa /4) 2 . By contrast, the four-band 
model includes the effect of the Coulomb interactions 
through its dependence on l. n 



In Table IV we show 



the ratios of the different Knight shifts and relaxation 
rates contributions. 



IV. THE STRONG INTERACTION LIMIT 



For strong interactions the four-band system in (2.1) 
is much more difficult to solve. Yet, it is still possible 
to highlight the qualitative features of the transferred 
hyperfine coupling interaction, specifically for the half- 
filled case. To obtain the strong interaction limit of this 
model we can perform the Gutzwiller projection 



H = P(H Q + H S + H N )P 



(4.1) 



which eliminates doubly occupied states in the Cu-3d or- 
bitals from the Fock space. The pro jection is effectively 
performed on all three terms of ( |4.l| ), which are treated 
on equal footing. As far as the first part PHqP is con- 
cerned, two possible superexchange processes are gener- 
ated, as shown in Fig. pi In the strong interaction limit 
(U a ^> \e v — e v >\,Ub ^> trfrj') the superexchange process 
in Fig. 0(a) is much more effective than the process in 
Fig. g(bj. For the basic system Ho, we only have to keep 
3 sta tes per unit cell, whereas for the four-band model 
(4.1), we end up with a system where we have to keep 27 
spin-degenerate local states per unit cell j with 4 tight- 
binding parameters t nv , for a half-filled system (exclud- 
ing doubly excited A, B-states, see Appendix |c|. For a 
doped system the number of states as well as the num- 
ber of possible transitions increases very fast, as has been 
shown for a tw o-ba nd model E3 A correct projection pro- 
cedure such as (|4.l|) becomes very difficult to handle, and 
one must resort to some approximations. In any event, 
in the vicinity of the half-filled case where the projec- 
tion can be explicitly used for the full Hamiltonian, we 
will analyze the differences between the predictions of the 



four-band model and those of the approximated Hamil- 
tonians. So let us restrict our analysis to the half-filled 
case where only virtual double occupancies of the copper 
site are allowed and where electron-nuclear interaction 
processes require that the initial and the final charge dis- 
tribution be the same. We deal with electron-nuclear 
interaction processes where effectively one local Cu-3d 
spin will be reversed and then relaxed by the thermody- 
namic fluctuations of the Heisenberg model. We decom- 
pose (Hq+Hs+Hn) into (L+K). L includes all local and 
K all kinetic con trib utions of the complete Hamiltonian 
H introduced in ( |2.l| ) . Then we can expand PHP on the 
basis of the unperturbed eigenstates of L and compute 
the projected local s-orbital spin operators like PSijjP. 
For the details we refer to Appendix |C] and discuss only 
the final results. 

First we analyze some relaxation processes for the oxy- 
gen atom. The process shown in Fig. m a) is a transferred 
(T) contribution proportional to 



F, 



B,T,(a) 



t B a 



.es - (e a + U a )_ 

For the process shown in Fig. 0(b), we include a part of 
the superexchange process to avoid double occupation of 
the copper site, and the contribution is proportional to 



B,T,{b) 



tabtBa 



X e a - eb)(es - Cb) 



(4.3) 



Then the lowest order contribution to the general form 
factor for the oxygen is given by 



F B (p) = 2F B , L + 2F B , T cos{pa/2) 



with 



Fbx = 

Fb,T = nB,T,(a)FB,T,(a) + n B,T,(b)F B ,T,(b) 



(4.4) 

(4.5) 
(4.6) 



projected Shastry 

n B,T,(i) denotes the combinatorial factor associated with 
all possible processes yielding a contribution F B t,(i)- 
The factor 2 for the left-right symmetry is not included 
in n BT u\. Like for the superexchange processes (Fig. @) 
some processes are forbidden due to the Pauli principle. 
However, since all energy levels are assumed to be spin- 
independent the related amplitudes F BT u\ are the same. 
For copper we also distinguish between the transferred 
(Fig. 0) and the local contributions (Fig. ||). The trans- 
ferred contributions are proportional to 



F 



A,T,(a) 



tAbtab 



(*A 



F 



A,T,(b) 



£b)(e J 4 + < 

tAbtab 



2e b - U b ) 

-\ 2 



F 



A,T,(c) 



A,T,(d) 



(e A - e b )(e A - e a - U a ) 

tAa 



e A ~ £ a - U a 

tAbtabtAa 

(^A ~ £b) 2 (ea - £6) 



(4.7) 

(4.8) 

(4.9) 
(4.10) 



whereas the local contributions are given by 



F 



A,L,{a) 



A,L,(b) 



A,L,(c) 



tAbtab 



(tA - £b)(tA + £a - 2e fc - Ub) 
2 



tAbtab 



(e A -e b )(eA - e a - U a ) 

tAbtab 



(eA - e b )(e A + e a - 2e fc )_ 
Then, the general form factor for copper reads 
F A (p) = 2F AtL + 2F AtT cos(pa) 

with 
F At L = n AtL ^ a )F A>L ^ a ) +n AL ( b )F A:L: ( o) 

+ n A,L,{c)FA,L,(c) + ■■■ 



F A ,T = nA,T,(a)F A ,T,(a) + nA,T,(b)F A ^T,(b) 

S v ' 

projected Mila— Rice 
+ ^A,T,(c)^A,T,(c) + n A,T,{d)FA.T.{d) + 

* v ' 

projected Shastry 



(4.11) 
(4.12) 
(4.13) 

(4.14) 
(4.15) 

(4.16) 



Let us now compare the predictions of the four-band 
model and those of the ID Mila-Rice or Shastry mod- 
els. 

Using the projected e xpres sion for the oxygen ins tead 
of the unprojected one (2.24), only process (b) in (4.6) 
contributes in t he strong interaction limit, whereas pro- 
cess (a) in (4.6) proposed by Shastry becomes negligible 



F, 



B,T,(a) 







(4.17) 



The form factor for the characteristic wave vectors (p = 
or p = tt/cl) is then reduced to 



Fb(O) = 2n B ,T,(b)FB,T,(b) 
F B (n/a) = . 



(4.18) 
(4.19) 



Since only the relaxation process of the oxygen nuclear 
spin contributes, which corresponds to p ~ , we re- 
cover the basic structure of the form factor of Shastry 
with modified amplitudes. Thus at half- filling, there is 
no fundamental diff eren ce for the oxygen between the 
general form factor (4.4) and the form factor proposed 
by Shastry fl2.27|). 



In the strong coupling limit at half filling, the follow- 
ing contributions to the form factor for copper are sup- 
pressed: 



F A,L,(b), FA,T,(b)> Fa,T,(c) 



0. 



(4.20) 



Thus the projected Shastry contribution (c) in (4.16 ) and 
one of the projected Mila-Rice contributions (b) in ( 4.16| ) 
as we ll as one of the projected local contributions (b) in 
(4.15) become negligible, and we end up with 



F a ,l = n A , Lj ( a )F A ,L,{a) + nA,L,(c)F A ,L,(c) + • • • (4-21) 
F A ,T = n Ai T,(a)F A ,T,(a) + n A ,T,(d)FA,T,(d) + ■■■ (4.22) 

for the local and for the transferred contributions to the 



general form factor (4.14), respectively. Thus the uniform 
part of the form factor is given by 

F A (0) = 4n At ( a )F Ay ( a ) + 27j j4 ,l,( c )-FXl,(c) 

+ 2n AfT{d) F A<T{d) + ... , (4.23) 

whereas the oscillating part reads 

F A (n/a) = 2n J 4,i,( c )-pA,i,(c) - 2n A ,T,{d)FA,T,{d) + ■■■ ■ 

(4.24) 

We used the fact that n A j J ^ a - ) — n At T,(a) — n A,(a) an d 
F A,L,(a) = F A,T,(a) = F A,(a)- The uniform part includes 
contributions which are absent in the ID version of the 
Mila-Rice and of the Shastry model. Furthermore, some 
terms proposed by Shastry turn out to be zero in the 
strongly interacting limit. For the oscillatory part the ef- 
fects are much more drastic. The transferred terms pro- 
posed by Shastry vanish in the strong coupling regime, 
whereas other transferred terms, which come from a com- 
bination of Mila-Rice and Shastry processes, contribute. 
Besides, the transferred terms proposed by Mila and Rice 
are canceled by the equivalent local terms. Hence, in 
ID, the general form factor differs both qualitatively and 
quantitatively from the form factors one would derive 
from the Mila-Rice or from the Shastry models. 



V. DISCUSSION AND PERSPECTIVES 

In this paper, we have analyzed the ID analogs of the 
hyperfine form factors proposed for NMR measurements 
of high-T c materials in the antiferromagnetic phase. We 
have focused on the situation where one deals with an 
antifcrromagnet generated by a superexchange process 
via an oxygen atom located at the midpoint between two 
copper atoms and where the Fermi contact interaction is 
one of the main contributions to the possible electron- 
nuclear interaction terms. We have investigated a ID 
Cu-0 model including four orbitals per unit cell, namely 
the Cu-3d and the 0-2p orbitals governing the ground 
state properties, as well as the Cu-4s and 0-3s orbitals 
describing the isotropic Fermi contact interaction. In ID, 
we were able to solve this model using only standard 
techniques without having to introduce any additional 
approximations for the hyperfine interaction term as pro- 
posed by Mila-Rice and by Shastry. Thus, we were able 
to compare our solutions of the four-band model with the 
predictions of the approximative models. 

In the low interaction limit, we have calculated the re- 
sulting temperature dependence of the Knight shifts K 
and of the relaxation rates 1/Ti for an undoped and for 
a doped system; in that limit the ground state is well 
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described by the strongly hybridized Cu-3d-0-2p anti- 
bonding band the width of which is large as compared 
to all Coulomb interactions. For both models, the four- 
band and the approximative one (Bulut model) , the tem- 
perature dependences are the same and show the typical 
power law behavior of one-dimensional interacting sys- 
tems (see Table |]J . Within this scope we have shown for 
the four-band model that for an undoped and a slightly 
doped system copper and oxygen behave completely dif- 
ferent for long distances or long times, when the tem- 
perature is low enough. The oxygen nuclei see only the 
Korringa-like contributions, since the antifcrromagnetic 
contributions are exponentially suppressed depending on 
the ratio of the characteristic length related to the charge 
gap and the doping. In contrast, the copper nuclei always 
see both contributions, the Korringa-like contribution as 
well as the antiferromagnetic one. This fundamental dif- 
ference between copper and oxygen vanishes gradually 
when the characteristic doping length or the character- 
istic thermal length becomes shorter than the length re- 
lated to the charge gap (the difference goes away abruptly 
when the system develops a gap in the spin sector). This 
solution is at variance with the prediction of the related 
approximate model, where for oxygen the antiferromag- 
netic contributions to 1/Ti increase with doping like 5 2 , 
whereas for copper they decrease proportionally to S 2 . 
Thus, the scenario where oxygen does not see the antifer- 
romagnetic fluctuations is realized much more effectively 
in the four-band model than in the ID version of the 
models proposed for the high-T c materials. In ID, such 
an unconventional scenario works, since even small in- 
teractions generate strong antiferromagnetic correlations 
due to the drastic reduction of the Fermi surface. 

We have also considered the strong interaction limit. 
Performing a Gutzwiller projection onto the four-band 
model without further approximations for the electron- 
nuclear interaction term, we computed the various pro- 
cesses contributing to NMR. Our analysis was limited to 
the insulating phase (Heisenberg model), since even in 
ID a full solution of the model for a doped system (t- 
J model with four orbitals per unit cell) is unavailable. 
In the strong interaction limit of the ID Cu-0 model, 
we were able to compare the form factors obtained for 
the four-band model with the predictions obtained for 
the approximate models (Mila-Rice model and Shastry 
model) investigating the different relaxation processes for 
the copper and oxygen nuclear spins. In this context, we 
have shown that neither the ID analog of the Mila-Rice 
model nor the ID analog of the Shastry model could de- 
scribe the strong interaction limit at half-filling. In con- 
trast to the usual assumption that only transferred con- 
tributions are relevant, we predict that both local and 
transferred contributions should be taken into account 
for describing the relaxation of the nuclear copper spin 
via an Cu-4s orbital. Furthermore, we have shown that 
for infinite local repulsions on the copper sites and small 
local repulsions on the oxygen sites, the contributions 
proposed by Mila-Rice and Shastry vanish. For the relax- 



ation of the nuclear oxygen spin we recover the basic idea 
of transferred hyperfine couplings with slightly modified 
amplitudes, but once again the contribution proposed by 
Shastry vanishes for infinite repulsion on the copper site. 

Both the strong and the weak coupling limits under- 
score the importance of keeping the full four-band model, 
at least in one dimension, in order to give an accurate de- 
scription of the NMR properties. The method we used 
in the present paper to tackle such a model can thus be 
extended in various directions. First, it can be applied to 
study specific models which have a structure similar to 
the model Cu-0 chain analyzed here. This is for example 
the case for ladder materials such as Sri4_ x Ca x Cu2404i. 
Analyses of the NMR material have so far been performed 
in terms of Mila-Rice-Shastry approximations. An anal- 
ysis retaining the full four-band model, with the spe- 
cific sympjjctrics of these ladder systems, is currently in 
progress.cll Other systems for which our analysis can be 
relevant are TMTSF and TMTTF alloys E At stoichio- 
metric composition they form an alternate stackEj Let 
us now comment on anisotropic contributions to K and 
to 1/Ti; these can be produced by a dipolar hyperfine 
coupling (see Appendix H) . They also stem from the spe- 
cific structural details of a given compound which may 
lead to an anisotropic form for the susceptibility: i n tha t 
situation the anisotropy of the p — component ( 3.37 ) 
will u sually be different from that of the p = 2fcp part 
( |3.38 ). In both the weak and the strong interaction lim- 
its, we find that - for l ow en ough temperature - 1/Ti is 
mostly determined by (3.3S), whereas K is proportional 
to (3.37). The experimental observation that 1/Ti is es- 
sentially isotropic and that K is anisotropic suggest that 
anisotropic effects are not too important for the p = 2kp 
contributions but do affect the p = terms. 

Another possible extension of our work concerns of 
course the two-dimensional systems. Although it is un- 
clear how much of the weak coupling approach remains 
valid in higher dimension, our strong coupling analysis 
can straightforwardly be applied to higher dimensional 
structures. The main difference in that case between the 
2D (or higher) and the ID study presented here comes 
from the symmetry of the various orbitals. In the case of 
a Cu-0 plane, in the presence o£,a Coulomb repulsion on 
the oxygen sites (Uo-2n ¥" )," th e rela ted amplitudes 
for the local processes (4.11) and (4.13) are not equal 
anymore, and a cancellation of these terms by symme- 
try arguments as assumed by M ila-Rice does not occur. 
Only the contribution like ( 4.12 ) will vanish by symme- 
try arguments. The transferred Mila-Rice contributions 
(4.7) via the 0-2p orbital, which always co st th e Coulomb 
energy Uo_2 P , and the local processes (4.11) have ex- 
actly the same combinatorial factor and the same am- 
plitude; thus the term (4.11) cancels out the term (4.7) 
for the antiferromagnetic wave vector. This suggests for 
^Cu-3d ~> °° that the antiferromagnetic contribution to 
the relaxation of the copper nuclei via an isotropic inter- 
action comes from local terms (see (4.13)) and fro m ne w 
transferred combined terms of third order (see (4.10)), 
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while the transf erre d cont ribu tions proposed up to now 
are absent (see (4.8) and (4.9)). 
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APPENDIX A: 



ANISOTROPIC HYPERFINE 
COUPLINGS 



Taking into account anisotropic hypcrfine couplings re- 
lated to the orbitals Cu-3d and 0-2p we have to replace 
(13) by 



m 



N 



js 



a 









c. 






b J j D bj 



(Al) 



The sum on ^ is over components of the diagonal hyper- 
fine tensors C^. 

In the weak interaction limit of the four-band model we 
can perform the same calculations as done b efore , and we 
will end up with the bosonized expression ( |3.13| ), where 
now the electron-nuclear interaction is given by 

h n = Y. fdx(c A \x Aa \ 2 + ci\x aa \ 2 )r(si + si a ) 

+ J2 J dx(C B \\ Ba \ 2 + C<\\ ba \ 2 )J<(S" a + S< a ) . 



(A2) 



Here we used the fact that S Ac 



and S Bc 






In general the explicit bosonized expressions for the spin 
part of S* and S^ a in ( |3.20 ) differ from S* a , but finally 
for a spin-symmetric model there will be no influence 
on the correlation functions. Thus only the coefficients 
are slightly modified and vary for the different directions 
<, = x,y, z. Formally the contributions to the ^-directions 
of the Knight shifts K^ u and Kq, as well as the contri- 
butions to the relaxation times T^ Cu and T£ are given 
by (3.35) and (3.36) performing the replacements 



|Aaq| 

|Asa| 



\^Aa\ + 7^-|Aaa| 

\^Ba\ + -pr~ \^ba\ 
^B 



(A3) 

(A4) 



in the expressions of the form factors defined in Table O. 

In the strong interaction limit of the four-band model 

the inclusion of anisotropic hyperfine interactions results 



0,3 



H'^ = J^(2C A F AtL + ClF a , L )r j Sl 

+ J^F AtT I](Sl J _ 1 + Sl J+1 ) 

+ £ ( C bFb,t + ClF btT )J^Sl } + Sl j+l ) 



(A5) 



js 



Here the new defined parameters which describe the ad- 
ditional couplings to the local Cu-3d spins are given by 
F a ,L = 1 an d F b}T = t 2 ab J_(e a — e b ) 2 , whereas all the others 
were defined in Section [V . For the copper atom the local 
contribution is modified, whereas for the oxygen atom it 
is the transferred one. 



APPENDIX B: THE SINE-GORDON MODEL 

At half- filling the spin part as well as the charge part 
of the Hamiltonian ( [3.13 ) are described by a sine-Gordon 
model H V SG = H% + H\j where 



m = 



2aU u 



(27ra) 2 



dxcoslVS^,, 



(Bl) 



For this model two different regimes exist depending on 
the value of the parameter K v . A massive regime {v m ) 
for 



2%u v {K v -l)< \U V \, 



(B2) 



where the perturbation of H]j is relevant, and a massless 
(r/ ) for 

2-ku v {K v - 1) > \U V \ , (B3) 

where the perturbation is irrelevant. 

1. Massive regime (i/ m ) 

When the cosine term is relevant, the conformal sym- 
metry is lost and the elementary excitations become mas- 
sive particles. To compute the correlation functions we 
can approximate the cosine term by 

m 2 r L 
H ™ = ^j o dx($„-($ u )) 2 , (B4) 

where the mass can be obtained from the exact solution 
of the sine-Gordon equation. For small U u one has 



f *K v \U v \ 



(B5) 



This Hamiltonian describes the fluctuations <5$„ of $„ 
about its mean value ($„) = 0. For such a system the 
Green's function (T T ^ lJ (r ly )^ u (0)) Um of the Laplace op- 
erator defined on the domain A v — [0 < u v r < u„/3, < 
x < L] is given by 
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G^(rv) = ^KoK(rv+a)] . (B6) 

_Ko is the Bessel function of zero order. 



2. Massless regime (i/ ) 

In this regime, the bare parameters are renormalizcd 
up to the fixed point values u v — > u*, K v — v if*, and 
t/j/ — » without changing the basic Luttinger Liquid be- 
havior of the unperturbed part Hq . For this model, the 
Green's functions for the unperturbed part regularized 
for large distances by R v and for short distances by the 
lattice constant a can be expressed as 



G v °(r u ) = ^ln[R u /( ri , + a)} 



or as the following limit 



G u °{r v ) = lim G"™{r v ) 



3. Correlation functions 



(B7) 



(B8) 



Typical spin-spin corre lation functions of the original 
fermions defined in (3.12) are combinations of exponen- 
tials of <&„. For a Gaussian model these functions can be 



expressed in terms of the Green's functions (B6) or (B7) 
depending on the phase vi, 

(exp[i7i$„(l)] . ..exp[i'y N $ v (N)]) I/i = 

>>B 7.^.(rH g -JE.TlC" 1 (C) | (B9) 



APPENDIX C: LOCAL STATES IN THE 
STRONG INTERACTION LIMIT 

The projected Hamiltonian (|j) is expressed as 

PHP = P{L + K)P , (CI) 

where L denotes the local system, whereas K includes all 
possible hopping terms of H. The eigenstates of L are 
given by 



;v 



|m,n 2 ,.. -,nj, ... ,n N ) = Yl 



(C2) 



j=i 



where n^ labels the local states n on site j. The local 
states and energies are shown in Table M. For simplicity 
we use the short notation 



|0i,0 



, u 2 , . . . , llj, . 



,mj,...0jv-i,0jv) = \rij,mi) (C3) 



(local ground state configurations are labeled by 1 j ) ) . 
The energy of such a state is given by 



E n , m = {N- 2)e + e n + , 



(C4) 



Now, we can expand the projection operator P onto 
the unperturbed eigenstates of L. Here for the half-filled 
case, we are only interested in the projection P onto the 
state |0) = n 7 =i |0j) with the energy Eq = Neo, thus we 
get 



« 



p = j:p 



where the first orders are given by 
P(°) = P P 

P {1) = PPQ + PqP 

P (2) = PPQQ + PQPQ + PQQP 

— (PpPQ 2 + PpQ 2 P + Pq 2 Pp) 



(C5) 



(C6) 



with 



Pp =Po 
Ppq =P KQ 



En-L 



Qo 



QP 



= Qz 



l 



Eo-L 



QoKP n 



Ppqq = PoKQ c 



1 



Eq-L 



Qo^Qo- 



1 



-Qo 



P 



QPQ 



P 



QQP 



= Qz 



= Qz 



l 



En-L 



E -L 
Q KP KQ —^-—Qo 



En-L 



En 



-QoKQo 



1 



En 



-QoKPo 



P P PQ 2 = PnKPnKQn 
Pp Q 2 P = PnKQi 



1 



(E - Lf 



1 



{En - L)- 



-QoKPo 



Pq 2 pp— Qo 



{En - L) 



T QoKP KPo 



(C7) 



The projection operator Q denotes 1 — Po- 

We can compute the projected electron-nuclear inter- 
action term PHpjP. The projection affects only the elec- 
tronic spins, and we have to evaluate projected local s 
orbital spin operators such as PS n jP. For example, the 
second order processes (see Fig. g(a) and 0(c)) are given 
by 



j(2) 



(P(°)+PW)S^(PW+P( 1 



(C8) 



Introducing (C§) and (C7) in M only 



Sg = PpqSa Pqp (C9) 

will contribute, because in a half-filled system the first 
hopping process brings the system out of the ground state 
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|0) and the second one brings it back to a possible ground 
state configuration |0). Here only the projections 



PPQ = 
Pqp = 



tAa 



Eq — E il3 

tAa 



Eo — E4 



>.:>, 



|0)(4„13 J±1 | 
|4„13 J±1 )(0| 



(CIO) 
(Cll) 



could generate finite matrix elements for second order 
contributions to the Cu-4s spin Saj- For the 0-3s spin 
Sbj it will be 



P P Q = 
P Q P = 



— tBa 

Eo — E12 

— tBa 
En — Eyi 



|0)(12 J | 
|12,)(0| 



Ppq = 
P Q P = 



tBa 



Eo - E 21Z 

tBa 



Eo — E' 



2,13 



|0)(2„13, +1 | 
|2„13 J+1 }(0|. 



Finally, the projected spins S r ■ are given by 



(2) 



m 



S A j — \\Aa\ Po (S aJ -i + S aJ + i) Pn 
^Bj = l-Vsal Pq (S a ,j + S aJ + i) Pq . 



with 



A/la = 
^Ba = ± 



tAa 



e A - £ a - U a 

tBa 

es - e a - U a 



(C12) 
(C13) 

(C14) 
(C15) 

(C16) 

(C17) 

(C18) 
(C19) 



Higher order contributions could be computed using the 
same procedure as for the above examples. 
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TABLE I. The lowest order temperature dependence of the 
Knight shifts and of the relaxation rates measured on the 
copper site (A) or the oxygen site (B) at half-filling (0) and 
for finite doping (5). Within a line, the undefined constants 
const, are the same for copper and oxygen. Temperature 
dependences are given up to logarithmic corrections. 





Copper 






Oxygen 




K° 

K n 


Fac x const. 
Fac x const. 






Fbc x const. 
Fbc x const. 




1/TS, 
1/Tr , 


(Fao.) 2 x const. 
(Fac,) 2 x const. 


x 


T 


(Fsa) 2 x const. 
{F Bo f x const. 


x T 




(Fac) 2 x const. 
(Fac) 2 x const. 


X 
X 


T 

rpK* 


(Fsa) 2 x const. 
(F Ba ) 2 x const. 


x T 
x T A '* 



TABLE II 
in the weak interaction limit at half-filling (7 = 0) and for small doping (7 



F va and F^ a are the non-oscillating and the oscillating contributions to the form factors of the four-band model 



8) for copper (A) and oxygen (B). The coefficients 



\ va come from the projection onto the lowest band (a) for small Coulomb interaction. Is and l mc denote the characteristic 
lengths related to the doping and the mass for a charge-gap system and K c is the Luttinger liquid parameter which controls 
the anomalous exponents of the correlation functions. Finally, Kq is the Bessel function of zero order. 



Copper 



Oxygen 






\\Aa\ 



2 I A., 



2|AAa| 4 cosh K c Ko (j 



(*)' 



|Asa 



2|Ab q | sinh 



K c Ko 



(*) 



TABLE III. The non-oscillating and oscillating form factors at half-filling (7 = 0) and for small doping (7 = S) for copper 
and oxygen. F^ 1 and FJ<f u denote the form factors obtained for the weak interaction Bul ut model with a Mila-Rice and 
Shastry like isotropic electron nuclear interaction term. The coefficients \ va are given by (3.4) and (3.5). By contrast, X V a is 



the characteristic coefficient of the Bulut model related to the overlap between the Cu-4s orbitals (A) the Cu-3d orbitals (a) 
or respectively between the 0-3s orbitals (B) and the Cu-3d orbitals (a) performed in real space. The other parameters were 
explained in Table ||. 



Copper 



Oxygen 



0Bu\2 



pBu 



Ol \Sui2i \ 2 
Z \ A Aa\ \Aaa\ 

s\x^\ 4 \x aa \ 4 (^-) 



8|A^| 4 |A aQ | 4 cos(7ra/Z 4 ) 2 cosh 



K c 



K c Ko 



(i) 



r >w Bu \ 2 w 

Z \ A Ba\ A ' 


(if;) ° 8|A|^| 4 |A aQ | 4 sin(7ra/2;,) 2 cosh 



KcKo (,£-)] (^) 
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TABLE IV. The ratios of the oscillating and non-oscil- 
lating contributions to the oxygen (rj — B) and the copper 
(j) — A) Knight shifts and relaxation rates at half-filling (0) 
and away from half-filling (S) calculated for the Bulut and 
the four-band model in the low interaction limit. The ratios 
are measured in units of their characteristic overlaps with the 
ground state (|A^*||A a| or |At,c<|). Here the oscillating con- 
tributions do depend on the characteristic length l mc and Is, 
whereas the non-oscillating contributions do not. 



Bulut model 



Four-band model 



( -0,<5 



l IT • 
I It 

1/T° 



TO" 



1 /T° 
'IB 

1/T<> 
> 1A 



VTf B 

TJrf A 



sin(7ra/2i,j) 2 
cos(7ra/l§)' i 



tanh K c Ko 



TABLE V. Possible local states per unit cell of the four 
band Hubbard model. The degeneracy is related to the pos- 
sible spin configurations for a local state labeled by |n,}. 



degeneracy 


local energies 




state 


2 


£7 = € A + 2e B 




|7.,) 


2 


eg = 2e A + e B 




|6,> 


1 


eg = 2<la + 2e B 




|5,> 


4 


ej = £A + 2e B + e a 




|4i) 


4 


eg = e A + 2e B + e b 




|3,> 


4 


e§ = 2eA + e B + e a 




|2.,> 


4 


ej = 2eA + £b + £& 




II7) 


2 


e = 2<la + 2e B + e a 




|o,-> 


2 


£l = 2e^ + 2e B + eb 




1.7 ) 


8 


e2 = e^ + 2e B + e a + e b 




|2,> 


2 


e 3 = e A + 2e s + 2e b + U b 




|3.,> 


8 


e 4 = 2ca + e B + e a + e b 




K,) 


2 


e 5 = 2e A + e B + 2e& + % 




|5,> 


4 


e 6 = 2ca + 2e B + e a + e 6 




|6,> 


4 


67 = e^ + 2e B + e a + 2ef, + [/& 


|7.,) 


1 


e 8 = 2e A + 2e B + 2e b + U b 




|8,> 


4 


eg = 2ca + e B + e a + 2e 6 + U b 


|9,> 


2 


em = 2ca + 2e B + e a + 2e& 


+ U b 


|io,> 


2 


en = ca + 2e B + 2e a + U a 




I11.7-) 


2 


ei2 = 2e A +e B + 2e a + U a 




|12.7> 


1 


eis = 2e A + 2e B + 2e a + U, 


I 


|13.,) 


4 


£14 = e^ + 2e B + 2e a + [/ a 


+ £b 


114.7) 


4 


eis = 2ca + e B + 2e a + [/ a 


+ £b 


IIS.,' ) 


2 


ei 6 = 2e A + 2e B + 2e a + [/< 


1 + £6 


II6.7 ) 


2 


ei7 = e A + 2e B + 2e a + U a 


+ 2e b + U b 


117.7) 


2 


eis = 2e A + e B + 2e a + U a 


+ 2e b + U b 


I18.7-) 


1 


619 = 2e A + 2e B + 2e a + [/< 


, + 2e b + U b 


|19i) 



1G 






J 


, 










i 






u b 


t 8b ~ 






LJ ^J 




8 

a 


r 


u. 




\ 


7 


e A 


r 


1 


■> 




\. 


J I 


J 




* 


A 


V 




r 


-\ r 


> 








\. 




7 



0-2p 
Cu-3d 

Cu-4s 
I 0-3s 



FIG. 1. Energy of the various orbitals. The Cu-3d and 
0-2p (resp. a and 6) orbitals are the relevant ones to de- 
scribe the electronic degrees of freedom. The Cu-4s and 0-3s 
(resp. A and B) should be kept to describe the coupling to 
the nuclear spin via a Fermi contact interaction. 




ja (j+l)a 

FIG. 2. Cu-3d and Cu-4s orbitals as well as the 0-2p and 
0-3s orbitals in the unit cells. The chosen signs of the wave 
functions determine the phase convention for the Hamiltonian 
and the signs of the various tight-binding parameters t vv i . 
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-X 

— r— 



9 I 



-X 

— h— 



1 I 



eJk) 



e.W 



FIG. 



H 1 1 1 1 1 1 

- -K o k F *- 

F F a 

3. Linearization of the spectrum close to the Fermi 
points. The momentum k is replaced by r Uf + q, where r — ± 
denotes the two possible directions. A > \q\ is an ultraviolet 
cutoff of the order of the bandwidth and the Fermi velocity is 

v F — d k e a (k)\k F - 



R 



r\a 



c -phase 



Q 




C 



a lm c /g r c 

FIG. 4. Crossover for the charge contribution to the 
spin-spin correlation function R^ a from the massive charge 
regime (c m -phase) to the massless charge regime (c -phase) 
at low temperature and for small doping. The behavior for 
copper (A) and oxygen (B) is quite different: For copper the 
amplitude remains finite for long distance, whereas the am- 
plitude for oxygen vanishes. 



(a) 



FIG. 5. Superexchange processes generating antiferromag- 
netic couplings between localized copper spins at half-filling. 
In contrast to the superexchange path (a) process (b) is sup- 
pressed, because it includes an intermediate state where it is 
necessary to pay the local Coulomb repulsion U a - The num- 
bers (l.,2.,3., and 4.) denote the sequence of the intermediate 
steps. 
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(a) 



(b) 



^^ 



-^4- .4h 



■Ore- 




0-2p 



Cu-3d 



0-3s 



FIG. 6. Possible spin flip processes of the nuclear oxygen 
spin via a Fermi contact interaction at half-filling, (a) is 
present in the ID version of the Shastry model, whereas (b) is 
another possible process which includes some superexchange 
contributions. At half-filling there are only transferred con- 
tributions (T), local processes (L) are absent. The numbers 
(l.,2.,3.,4. and 5.) denote the sequence of the intermediate 
steps. 



(a) 



(c) 




FIG. 7. Possible transferred spin flip processes of the nu- 
clear copper spin like PS\-P via a Fermi contact interaction 
at half-filling, (a) and (b) appear in the ID Mila-Rice model, 
(c) pertains to the ID the Shastry model , whereas (d) is a 
combination of both. 



(a) 



£,.- 



e __ 



e„- 




(c) 




^^0-2p 



Cu-3d 



Cu-4s 



3. 3. 

FIG. 8. Possible local spin flip processes of the nuclear copper spin like PS^jP via a Fermi contact interaction at half-filling. 
All processes include the 0-2p orbital as an intermediate state and are of the same order as the transferred hyperfine coupling 
processes proposed by Mila and Rice. 
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